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Abstract 

In this paper the hybrid-NLIE approach of [38] is extended to the ground 
state of a D-brane anti-D-brane system in AdS/CFT. The hybrid-NLIE equa¬ 
tions presented in the paper are finite component alternatives of the previously 
proposed TBA equations and they admit an appropriate framework for the 
numerical investigation of the ground state of the problem. Straightforward 
numerical iterative methods fail to converge, thus new numerical methods are 
worked out to solve the equations. Our numerical data confirm the previous 
TBA data. In view of the numerical results the mysterious L = 1 case is also 
commented in the paper. 


1 Introduction 


In this paper in the context of AdS/CFT [U [2, 3] we study numerically the ground 
state of a pair of open strings stretching between two coincident _D3-branes with 
opposite orientations in S 5 of AdS$ x S' 5 . The main motivation for the study is that 
according to string-theory the ground state of such a configuration is expected to 
be tachyonic for large values of the’t Hooft coupling [4]. In our work we rely on the 
perturbatively discovered and later ’’all loop conjectured” integrability [5] of both the 
AdS$ x S' 5 super-string and the dual large N gauge theory. For string configurations 
with D-branes integrability enabled one to describe string configurations ending 
on different types of D-branes as 1-dimensional integrable scattering theories with 
boundaries P0EIEI ITU] , This formulation of the problem makes it possible to 
go beyond the approaches of perturbative gauge and string theories being valid for 
small and large values of the ’t Hooft coupling respectively, and to determine the 
exact spectrum of the model at any value of the coupling constant. However, even 
with the help of the powerful techniques offered by integrability, the exact analytical 
solution of the problem is not possible. Remarkeble analytical results are available 
in the small mmmmmm and large [116., 1?!, [T9J IB] coupling regimes, but 
the determination of the spectrum at any value of the coupling constant can only 
be carried out by high precision numerical solution mmm of the corresponding 
nonlinear integral equations. 

In our paper we consider the case, when the two D3-branes are giant gravitons 
[23] , namely they carry N units of angular momenta in S' 5 . If the S' 5 of AdS$ x 
S' 5 is parametrized by three complex coordinates X, Y, Z satisfying the constraint: 
|A"| 2 + |Y| 2 + |Z| 2 = 1, then our D3-brane and anti-D3-brane are given by the 
conditions Y = 0 and Y — 0 respectively. They wrap the same S' 3 , but with 
opposite orientation. As a consequence of Gauss law such a system can support 
only even number of open strings. For this reason we study the minimal number 
of allowed open strings, a single pair, ending on our D-brane anti-D-brane ( DD ) 
system with open string angular momenta L and L'. 

On the large N gauge theory side a Y — 0 brane is represented by a determinant 
operator 121 composed of N copies of the field Y: 

Oy = det Y = y„»‘ ■ • • Yl« (1.1) 
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where a, and 6* are color indices and e is a product of two regular epsilon tensors 
= e ai '" aN eb 1 ..-b N ■ The local operator corresponding to an open string ending 
on a Y = 0 giant graviton can be obtained from (11.11) by replacing one Y held with 
an adjoint valued operator W [ 25] : 

Or = C-CU •'' Y% : ; . (1.2) 

The gauge theory description of a pair of open strings stretching between two D- 
branes is given by a double determinant operator, such that the string insertions W 
and V connect the two determinants of the Y heldqj]: 


n w,v 

U Y,Y 


ai—aAryfti _ _ _ y6 jv-i ci -cn ydi _ _ _ y d N _ 1 vyijv y£>jv 

1 a i 1 a N -1 ' z di---dpf 1 ci 1 * * Civ_i ajv ^ cjv 


(1.3) 


Unfortunately, the precise gauge theory dual of the Z5Z)-system of our interest is 
not known. In [4] it was approximated by a double determinant operator similar to 
(11.31) . but in one of the determinants the Y fields are replaced with Y heldj^j: 


n w,v 

U YY 


ai-.-aivy6i _ _ _ ybjv-i ci---c N yd\ _ _ . yrfjv-i VW^iv V? 6 iv 
t 6i..-bjv 1 a i 1 a N -1 ^di-'-d^r 1 ci 1 cn -i y v y cn ' 


(1.4) 


For the ground state the insertions are W = and V = Z L/ respectively. Based 
on one-loop results the planar dilatation operator is expected to act independently 
on the two words W, V corresponding to the open string states jUj : 


— Ab ar e [(9^ V ] + 5A[Wi'y] + 5A[Vyy]. (1-5) 


This observation allows us to apply the boundary Thermodynamic Bethe Ansatz 
technique (BTBA) [25] to each open string separately. The necessary ingredients of 
this technique are the boundary reflection factors [SJ EH EH E9] and the asymptotic 
Bethe equations of the problem [4], Unfortunately, apart from some very special 
cases [30, 31], [32], it is still unknown how to derive BTBA equations for a general 
non-diagonal scattering theory in the context of the thermodynamical considerations 
of [25] • This is why in [4] the Y -system [33, M105J and the related discontinuity [36J 

1 The ground state of such string states is BPS. 

2 According to the argument of [3] the correct state might have other structures involving the 

fields Y and Y, but should be similar to the double determinant form m and the mixing with 

other fields seem to be suppressed at large N 
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equations supplemented by analyticity assumptions compatible with the asymptotic 
solution BZJ were used to derive BTBA equations for the nonperturbative study of 
the ground state of the DD -system [4j. 

The BTBA description of the system is an infinite set of nonlinear integral equa¬ 
tions. The numerical solution of the equations [1] showed that the ground state 
energy is a monotonously decreasing function of the coupling constanto g. The ana¬ 
lytical investigation of the large rapidity and large index behavior of the Y -functions 
of the BTBA revealed that the usual BTBA description of the system breaks down 
when the energy of an open string state with angular momentum L gets close to 
the critical value: E C (L ) = 1 — L. This point was interpreted in [0] as a transition 
point where the ground state becomes tachyonic. Approaching the critical point 
the contribution of infinitely many Y -functions must be taken into account to get 
accurate numerical result for the energjQ. This fact suggests reformulating the finite 
size problem in terms of finite number of unknown functions. The possible candi¬ 
dates could be the FiNLIE [06], the quantum spectral curve (QSC) [47» M] or the 
hybrid-NLIE (HNLIE) [38]J formulation of the problem. Since at present it is not 
known (not even for the Konishi problem) how to use the analytically very efficient 
mm QSC method for numerical purposes, we choose the HNLIE method to re¬ 
formulate the finite size problem of the _DZ)-system. In this paper we transformed 
the infinite set of boundary TBA equations [4] into a finite set of hybrid-NLIE type 
of nonlinear integral equations. We perform the extensive numerical study of these 
type of equations in order to get as close to the special Ebtba = 1 — L critical point 
as it is possible. 

Our numerical results reproduce the numerical evaluation of the boundary Liischer 
formula [271 [59 ] in the linear approximation, and the numerical BTBA results of |T] 
as well. These numerical comparisons give further numerical checks on the hybrid- 
NLIE technique of (38J. Unfortunately, as g increases new local singularities enter 
the HNLIE formulation of the problem. Thus we could not approach very close to 
the critical point. Nevertheless, in the range of g where physically acceptable numer- 

3 Throughout the paper the relation between g and the ’t Hooft coupling A is given by: A = 
4ttV. 

4 This means that the usual truncation procedure for solving the infinite set of TBA equations 
is not applicable to such a system 
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ical results were obtained, the HNLIE results could give higher numerical precision 
than that of the BTBA and also some interesting facts conld be read off from our 
numerical data. 

During the numerical solution of the HNLIE equations straightforward numerical 
iterative methods failed to converge, thus new numerical methods were worked out 
to solve the equations. 

The ground state of the L — 1 state is a very special case, since there the critical 
point is right at g — 0 and so far neither perturbative field theory computations 
nor the boundary Liischer formula could provide a finite quantitative answer to the 
anomalous dimension of this state. On the integrability side the HNLIE approach 
allows us to get some numerical insight into this problem. 

The outline of the paper is as follows: Section 2. contains the HNLIE equations. 
In section 3. the numerical method is described. In section 4. the numerical results 
and their interpretation is presented. Section 5. contains some comments on the 
mysterious L — 1 case and finally our conclusion is given in section 6. Various 
notations, kernels of the integral equations together with the necessary asymptotic 
solutions are placed in the appendices of the paper. 

2 The HNLIE equations 

In this section we transform the previously proposed BTBA equations of [4] for 
the ground state of our D-brane anti-D-brane system to finite component hybrid- 
NLIE equations. For presentational purposes we group the equations into 3 types. 
There are TBA-type equations, horizontal 577(2) hybrid-NLIE type equations, and 
vertical SU (4) hybrid-NLIE type equations. They together form a closed set of non¬ 
linear integral-equations, which are solved numerically in this paper. As it is usual, 
structurally the equations consist of source terms plus convolutions containing cou¬ 
pling dependent kernels and nonlinear combinations of the unknown functions. The 
objects appearing in the arguments of the source functions are subjected to quan¬ 
tization conditions, but similarly to the boundary TBA description [3], due to the 
u —> —u symmetry of the problem they are tied to the origin of the complex plane, 
thus extra quantization conditions are unnecessary to be imposed, since they are 
automatically satisfied by symmetry. Since these source term objects have fixed 
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positions their positions are exactly the same as that of their asymptotic counter¬ 
parts. This fact saves us from the tedious computation of the source terms, since if 
we take the difference of the exact equations and their asymptotic counterparts the 
source terms cancel from the equations. To be pragmatic and save time and space, 
the equations will be presented in such a difference form. Thus for any combina¬ 
tion / of the unknown functions, we introduce the notation 5f{u) = f(u ) — f°(u), 
where f°(u) is the asymptotic counterpart of /. Having introduced this notation, 
we start the presentation of the equations by the TBA-type part. For the labeling 
of the Y -functions we use the string-hypothesis [30] based notations of [31] . For the 
presentation of the equations a few more notations need to be introduced: 


-2± = log 


r 


1 


= log 


1 + 


'^Tn\vw J 

For later numerical purposes we re-parametrize log Yq by the formula: 


t(u) = tanh(^p). 


( 2 . 1 ) 


log Yq(u) = -2 L log ^_ Q ’ +\ogy Q (u)+CQ+E log ( u 2 + U__-J , Q = 1,2,... 

3 ( 2 . 2 ) 

such that cq is the constant value of log Yq at infinity and £ is minus twice the 


energy^: £ = —2 E 


btba ■ From the TBA equations of the problem [3], it follows that 
5cq = cq — Cq = 5c is Q- independent, and for small g , log i/Q is a smooth deformation 
of its asymptotic counterpart, such that 5 log yQ tends to zero at infinity^ 

Using this decomposition the following notations are need to be introduced: 


Lq = log(l + Y q ), 5 R q = log(l + Yq) - 5\ogy Q . 


(2.3) 


Then the TBA-type equations take the form: 

<51og Y m \ vw 5 log [(1 T Uti+i|i«u) (1 T 5 m — i|i«u)] lo§(lT Y m +1 jics , 2 Y m ^ p g 2, 

(2.4) 


5 The log multiplier of £ in (12.21) is chosen not to modify the constant term in the large u behavior 
and to satisfy = - Sq S - q — 5 6 ®o+ > which is the LHS of an important Y- 

J l'Q-l'Q+1 Y° + 1° W-lSQ+1 Vq Vq 

system equation divided by its asymptotic counterpart. 

6 log Yq cannot be considered as smooth deformation of log Yq, because log Yq —log Yq ~ e log |u| 
diverges for large u at any g. On the other hand log yQ — log yQ is small for any u at small g and 
tends to zero at infinity. 
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5 log Y ,\ vw = 5 log (l + Y 2 \ vw ) *s - log(l + Y 2 ) k s + 5 log 


1 - Y_ 
1 -Y + 


(fiog Vq — 2 * s — (8 Rq_i + 5Rq + i * s), Q > 2 , 

Y pq—2 

siogY = -J]iog(i +Y Q )* K Q ,-a(K Q ,). 


* s, (2.5) 

( 2 . 6 ) 

(2.7) 


Q =1 


S\og(Y + Y_) = 2hlog 


1 + Y\ | t 


P0-2 


*» +£ log(l + Y q ) * [-AT, + 2/i'O , 1 *«] 
Q=1 


L1 + Y 1{ 

- n(K Q ) + 2(l(K%*s) 

For Y\ the modified hybrid form m of the BTBA equations is used, 


( 2 , 8 ) 


8 log yi =2<5 log(l + Yi | vw) S ~k Kyi 
1 - Y_ 


— 2<5 log 


PO-2 


i-n 


-k S -k 


/o'!, + 2^if_ * K y } + 2M? + * Kf 


(2,9) 


+ V log(l + Vq) * [a«‘ 2) + 2s * A'®" 1 ' i 

0=1 

+ !J(A« 1 2| ) + 2!J( S *A«- 1 ’ 1 ), 

where p 0 is the index limit starting from which the upper part of the TBA equations 
is replaced by an At/ (4) NLIE of [38] (See figure CD). For any kernel vector appearing 

OO 

in the TBA equations h2(/C q) denotes the residual sum ]T) Lq*/Cq, and following 

Q=po-l 

the method of 021 for po > 4 it can be expressed by next to nearest neighbor Y- 
functions as follows: 

m Q ) 5 Rp Q —i k (J l k KjpQ—2 $RpQ—2 k (7^ k KLp^—\ 

T 25r pq— 2 ^ ^ KpQ—2 2 5rp 0 —3 ^ si ^ (2.10) 

where r m = log(l + Y m \ vw ), the kernels s, s 1 , ai are hyperbolic functions 

ivgu 


s(u) = 


9 


4 cosh ^ ’ 


/ X 1 ,u. / x 9 cosh' 

1 W= 2 S( 2 l = 


( 2 . 11 ) 


while the other TBA kernels can be found in appendix A. As a consequence of the 
re-parametrization (I2.2|i the two constants 5c and e also become part of the set of 
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( 2 . 12 ) 


equation^ 




2vr 


Q =1 


du 27t du 


5c =2<51og(l + Yi\ vw ) * s * CK y i 

"i -y_ 


2d log 


PO-2 


1-K, 


-k S k 


CK “ + 2i25?_ * CAK 1 + 2<S_S < ’ + * CA'f 


(2.13) 


+ V log(l + y„) * [c/f®, + 2s * CAf„‘. 1 

Q=l 

+ SlfCAT®,) + 211(s*CA'0- 11 ), 

where for any kernel 1C: ClC{u) denotes the constant term in the large v expansion 
of /C(n, v)u As we mentioned —e/2 is the TBA energy, thus (12.1 2j) gives the energy 
formula in our formulation of the finite size problem. The asymptotic forms of the 
Y -functions necessary for the formulation of (I2.4H2.13^1 are listed in appendix D. To 
close the discussion of the TBA-type equations we note that equations (12. 7p and 
(12.8(1 determine Y± up to an overall sign factor. The sign factor can be fixed from 
the asymptotic solution and its value is —1. Thus the fermionic Y-functions can be 
expressed in terms of the LHS of (12.7(1 and (12.8)1 by the formula: 


Y* = 


;V+Y_±i: 


Y_ 

Yl 


(2.14) 


The horizontal SU(2) wing of the TBA is resumed by an S77(2)-type NLIE 
which in our case takes the form: 

51og(-6) = si 1 ” 7 8 9 ! *5log(l + Yl,™) + G*S log(-l - b) - G [ " 271 *5log(-l - b), (2.15) 
5log(-6) = * 5 log(l + Y llw ) + G*5 log(-l - b) - G™ *5 log(-l - b), (2.16) 


hlog Yxi^ = 1 ^hlog(—1 — b) + ^ *<51og( —1 — b) + <51og 


l-i 


is, (2.17) 


where 0 < 7 <l/ 2 isa contour shift parameter, the kerne 
the asymptotic solution for b and b is given in appendix D@ 


G is given by (IB.611 and 
The upper SU{ 4) NLIE 


7 Here the * notation means simply integration from — oo to oo. 

8 Here we note that only the dressing kernel has logarithmically divergent term in its large v 
expansion, all the other kernels has either constant term or they simply vanish at infinity. 

9 In practice b and b are complex conjugate of each other. 
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of [38] is attached to the TBA equations at the p 0 _ hli node. The upper NLIE is for 
12 complex unknown functions: 6 a and cIa, A = 1,6. They are combinations of 
the T-functions of the upper wing 577(4) Backhand-hierarchy [3S] - Their relations 
to the unknowns introduced in [38] are given by (IB.3IIB.jp in appendix B and their 
asymptotic forms are given in appendix C. Using the notation Ba = 1 + 6a and 
Da = 1 + (Ia , the equations they satisfy take the form: 


<51og6 a = ^(GVb)aa' * 61og Ba' + ^(G7d)aa' * 61og Da> + Ea, 


(2.18) 


A' 


A' 


(hog d A = ’y^ J {GdB)AA' k tflog B A ' + Y J (. G dD)AA'*8\ogD A ' +E a , (2.19) 

A' A' 

where the kernels are given in (1B.5HB.12p . The shifts in the kernels which is equiva¬ 
lent to fixing the lines on which the NLIE variables live, are chosen in a symmetrical 
way and fixed as follows: 



i(—9,-1,5,—3,3,1), 
^(-5,1,9,-3, 3,-1). 


( 2 . 20 ) 

( 2 . 21 ) 


This choice satisfies the constraint inequalities of [38J and satisfy the relation 
7 = —Wit] with WX given by (1C.211) . Its advantage is that choosing the (7 = 0 
asymptotic solution from appendix C to formulate the equations, the 6- and d-type 
variables are related in a simple manner: 


b(-u) = md(u). (2.22) 

In practice this reduces to half the number of 5(7(4) NLIE variables. The vectors 
Ea and Ea are conjugate to each other and they give the TBA input into the upper 
NLIE. To give their form we introduce the notations: 


m = Y [ei] 

'1 pO— l\vw^ 

Vi ~ yJSi| w , £ 1 ~ ~ e 3 - 

7 

_ 12 ’ 

(2.23) 

*<$log(l + 771 ), 

E 3 = * 61og(l + fji), 

e 5 

= E e = 0, 




(2.24) 

E 2 = -S 1 k dlog(l + T] 1 ) - -s 1 3 ] k dlog(l + r/i) + £2 + i y> 2 , 

(2.25) 








1 1 1 

E 4 = - -s^ 1 * 51og(l + 771 ) - -S [ f ] ~k <flog(l + fji) + s 4 + i ip 4 , 


where 


£2 (u) = — / dv SRp 0 (v) 


£±{u) = 27r / dV SRpo ^ 


1 | 1 
» - , ; - W, M V ~ ife J ’ 


1 1 I 

n-v-i- a + u + v-i-J’ 


(2.26) 


(2.27) 


(2.28) 


<P 2 (u) = [ dv hlog(l + Vi(v)) ( ip{u -v + ^-) + <p(u + V - |^) 


(2.29) 


p 4 (u) = [ dv <flog(l + r)i (v)) { <p(u -v + ^-) + p(u + u - |^) 


(2.30) 


( \ 9 f.,,1 iug . ,,1 , iug ^ irgu 

vt u ) = W(4 “ ~r) + —) - 71 tanh ( 2 ) 


(2.31) 


The last set of equations gives, how the upper NLIE variables couple to the TBA 
part of the equations. 


5\ogYp 0 _ 2 \ vw = s [ ~k hlog(l T Tj\) + s*hlog(l + Yp 0 — 3 \vw) ~ s*L P0 _ i, (2.32) 
hlog T]i = <s [ -1+ei- T' 1 ] * ^l Q g _|_ s [l+ei-m] * ^l 0 g £) 1 _ s [ e i- 72 ] * ^J Q g _g 2 
+ s [ei] * hlog(l + Y p0 - 2 \ vw ) — * L po , 

Slog fj! = s [_1+e3_731 * flog B 3 + s [1+e3 ~ m] * flog D 3 - s [e3_,?2] * hlog D 2 

(2.34) 

+ s N * flog(l + Y P0 _ 2 \ vw ) - * L Po , 

flog y P0 = s 1 - 1 -^ * [Slog b 2 - flog B 2 \ + s [1 -' ?2] * [Slogd 2 - flog D 2 ] 


+s [ " 73] * flog ^ + s [ -' n] * hlog ^ + a W ★ flog 

6 3 di r] 1 

+S [_e3] * flog -r — - s* 5 R„ 0 _ 1 , 
hi 

where 62 and c ?2 are from the re-parametrization of b -2 and d 2 - 


h(u) = V — 


xi- po+72l (u) 


(2.34) 


(2.35) 


/ [i f , • 72 ^Po - 1 \ . . tt] (fcl - 
exp | £ loglw + 2 --- ) +7 2 + yj 

(2.36) 
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2 L 


d 2 {u) = g 


x [ r +V 2 ] (u) 


exp 


log iu + i 


. m + Po + 1 
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. vr 
1 2 


+ - 77 - } d 2 (u), 


(2.37) 

with g = ±1 being a global sign factor. Similarly to the definition of t/q, also 
here the benefit of using b 2 and d 2 is that, for small g, log & 2 and logd 2 are smooth 
deformations of their asymptotic counterparts, and in addition hlog & 2 and hlogd 2 
vanishes at infinity, which is necessary for the convergence of certain integrals. The 
decompositions fl2.36p . fl2.37p are chosen to be compatible with the functional relation 
& 2 -72 ^ d 2 ~ 7 '^ = Y p0 in ref. [38]. Equations (I2.4p - (l2.35p constitute our complete set of 
nonlinear integral equations, which governs the finite size dependence of the vacuum 
of our D-brane anti-D-brane system. 


3 The numerical method 

ffere we describe our numerical method for solving the hybrid-NLIE equations pre¬ 
sented in the previous section. During the iterative numerical solution of the equa¬ 
tions we faced with very serious convergence problems, which forced us to work out 
such a method that overcomes all the difficulties emerged. Our numerical method 
can be applied to solve other type of nonlinear integral equations as well. The power 
of the method is shown by the fact that numerical convergence was reached even in 
such cases, when the solution was physically unacceptable. 

The numerical method consist of two main steps, namely: 

• Discretization of the equations 

• Iterative solution. 

The first step involves the discretization of the unknown functions and kernels, 
furthermore the discrete approximate representation of the convolutions. Having 
carried out the appropriate discretization of the problem, the equations are consid¬ 
ered as large nonlinear algebraic set of equations. Thus eventually instead of integral 
equations we solve discrete algebraic equations. In this paper we will present two 
methods to solve them numerically. 
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3.1 Discretization of the problem 


The discretization serves two goals. First it allows us to reduce the numerical prob¬ 
lem from solving integral equations to solving algebraic equations. Second choosing 
the discretization points appropriately it reduces the number of degrees of freedom 
as much as it is possible to reach the desired numerical accuracy. In our actual 
numerical computation instead of u of section 2. we used the new rapidity u —> 
because with such a scaling almost all the rapidity difference dependent kernels 
become g independent. Thus for example Y±{u) will be defined in [—2g,2g\. To 
decrease the number of discretization points the u —> — u symmetry of the problem 
is exploited. This means that the E-functions are to be discretized only on [0, oo] 
or [0, 2 g\ and as for the NLIE variables it is enough to discretize the b- and d- type 
variables on [0, oo]. Since we do not want to introduce any cutoff in the rapidity 
space first we transform the u £ [0, oo] interval to a finite interval t £ [0,5(a)] 
through the transformation formula: 


u{t) 



B(a) = 2- + 1. 
a 


(3.1) 


This formula is chosen such that the branch point 2 g corresponds to t — 1 for any 
choice of the parameter a, where a is a global scaling factor which changes from 
unknown to unknown. We chose the values as follows: for Y\\ w a = 1, for b and b 
a = 2, for Yq and Yq_ i\ vw a = Q, for r/i and fj i a = p 0 , and finally for the b- and 
d-type NLIE functions a = po. These values are chosen to preserve the smooth- 
nesa 10 l of the transformed functions in the finite interval. After this transformation 
all of our unknown functions live on a finite interval. To discretize them we used 
piecewise Chebyshev approximation. This means that we divide the finite interval 
into subintervals and on each subinterval the functions are approximated by a given 
order Chebyshev series. The choice of subintervals is not equidistant. The subinter¬ 
vals are placed more densely around the branch points, since the function x(u/g), 
which governs the decay of the massive Yg-functions, has the largest change around 
this point. The advantage of the Chebyshev approximation is that if the function 
is smooth enough on the subinterval, the coefficients of the Chebyshev series decay 


10 In our terms the lack of smoothness would not mean discontinuity, but the presence of rapidly 
changing parts and peaks. 
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rapidly and the order of magnitude of the last coefficient allows us to estimate the 
numerical errors of the procedure. Now we describe the discretization method in 
more detail. Our functions are defined on either [0, B(Q)\ or on [0,2 g]. This is why 
two type of subinterval vectors are defined Aq and A±, such that the endpoints of 
the subintervals of [0,5(Q)] are put into the vector Aq and the endpoints of the 
subintervals of [0,2g] define A±. Let l k be the order of the Chebyshev approxima¬ 
tion, then using the general rules of the Chebyshev approximation, a given function 
f{t) is approximated in the kth subinterval [A k -i,Ak] as: 


II. 

m - E 

3 = 1 


(fc) rri 

c ) T, 


3 -1 


ft — \{Ak + 

V \{A k — Ak-i) J ’ 


t e [A k _i, A k \, 


(3.2) 


where now the vector A stand for either Aq or A±, furthermore Tj_i are a slightly 
modified Chebyshev polynomial 


TM) = 


Tj(u ) if j > 1, 

§ if j — 0, 

with Tj(u ) being the jth Chebyshev polynomials. The coefficients are the 
Chebyshev coefficients of the function /, which can be computed from the sampling 
points of the Chebyshev approximation: 


J 2 
by the simple formula: 


t\ ^ — -^{A k — A k _i) C W(Z fc ) + — (A k + A k _i), i — 1, ,.,/fc 


(3.3) 


Jfc) = ± 

j h 




j 0 = 1 


where c^\l k ) are the zeros of the l k order Chebyshev polynomial: 


C«(4) = — COS 

and C k l is given by: 

C k)i = cos 


71 

— \ l — ~ 


Ti k (c (i \lk)) = 0, i = 1, ...,h 


7T 


r (i-l) 


i,k e {!,...,4}- 


(3.4) 


(3.5) 


(3,6) 


11 This slight modification is only to write the approximation series m in a more compact way. 
12 The Chebyshev polinomials are defined by the formula: Tj{u) = cos(j arccosit), j = 0,1,2... 
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In our method the next step is to formulate the convolutions and the equations 
themselves in terms of the discrete values of our functions. Here will sketch the basic 
idea in some typical scenarios appearing in our equations. Then its application to 
the concrete unknowns and kernels of the problem is straightforward. If one takes 
the equations at the required discretized points t) the following typical pattern 


arises: 


F(u(t { p)) ~ / dv' L{v') K s (v', + ..., 


d k 'h 


(3.7) 


where Ii s (u, v) = K(u, v ) + K(—u, v ) is the symmetrized kernel to exploit left-right 
symmetry of the problem for reducing to half the number of variables. F(u(t^ ^)) is 
intended to modelize the variables in the left-hand side of the equations taken at the 
discretized points of the transformed variable t and L(u ) s tands for some nonlinear 
combination of some unknown function of the equations^. If L(u(t )) is discretized 
by a subinterval vector A of [0, B(a)\ , then the numerical approximation of the right 
hand side goes as follows; 

• First the integration variable is changed from v' to t, 

• then on each subinterval L(u(t )) is approximated by its Chebyshev series, 

• finally the integration is carried out and the convolution is expressed in terms 
of the discretized values of L(u(t)). 

The final approximation formula takes the form: 


£(A) i k 


f dv'L(v')K 3 (v',u(tf } ))~Y i Y i L kJ (££$ 

n k=l j= 1 V k jo=l 


irk, j 0 
h-J+lj'o IK 'k',j' 


(3.8) 


where = L(u(t ( J k ' > )), C(A) denotes the dimension of A and F k k A-, is the discretized 
convolution matrix given by the formula: 


Ak 


dt 


K k ^, = aB{a) / 


A k ~ 


t — b(A k + Ak-i) 


\{A k 


Ak- 


k-1 


*■*(«(*), 


(3.9) 


13 For example in the TBA-part F(u) can be thought of as log Y(u) and L(u) can be log(l+T(u)) 
for any type of Y. 
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In this manner a convolution is reduced to a discrete matrix-vector multiplication. 

The other type of typical convolution is when the integration is taken from zero 
to 2 g. In certain cases the function L(u) has square root behavior close to the branch 
pointj^l. For such functions the truncated Chebyshev series does not give accurate 
approximation. In these cases not the function L{u) is approximated, but that part 
of it which remains after the elimination of the square root behavior. Namely, we 
write L(u) = \J ^g 2 — u 2 L(u), then L{u ) is approximated by a truncated Chebyshev 
series and finally the approximate discretized form of the corresponding convolution 
is very similar to (j3.8j) : 

c{A±) i k / i k \ 

dv L(v) K s (v, U { P) T E V& , (3.10) 

k =1 j=l \ k jo=l J 

where Lkj = L(v'P) and /Cjjh, is the square root factor modified version of (13.91) : 

A± : k 

= J dv yj 4 g 2 - v 2 fj_ 1 

A±,k-i 

Here depending on the left hand side of the equation can stand for u(ty' > ), 
t G [0, B(a)} for some a, or it can denote the sampling points on [0, 2 g). 

Applying our discretization technique to all unknowns and convolutions of our 
equations, we can reduce the integral equations to a discrete set of nonlinear alge¬ 
braic equations. However, the transformation from integral equations to algebraic 
equations is obviously not exact. The typical error comes from the fact that on each 
subinterval the Chebyshev series is truncated, so the magnitude of the typical errors 
in our numerical method is governed by the neglected terms of the Chebyshev series, 
which can be approximated by the magnitude of the last Chebysev coefficient. In 
our case this is typically somewhere between 10 -5 and 10~ 6 . 

The last step of our numerical method is the iterative solution starting from the 
asymptotic solution. 

^ ' < ^ -j^ y J-— y — 

14 Such typical combinations are log 1 _ y ~ and log 1 f . 


V — 1) 


K s {v,u [ p). (3.11) 


/ 
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3.2 The iterative solution 


Here we will describe two methods to solve our integral equations iteratively. Since 
our actual equations have very complicated form, we will describe our methods using 
a model example, which has similar structure to our equations. 

Let the model equations take the fornJ^I: 

log |/a = fa + Gab*\og(l + y b ), (3.12) 

where G a b are some kernel matrices, f a are some source terms and y a s are the 
unknown functions of the problem. The solution of (13. 12j) is expanded around the 
asymptotic solution and the equations are formulated in terms of the corrections. 
To fix the conventions, the correction functions 5y a are defined by: 

ya = ya( 1 + S Va)- ( 3 - 13 ) 


As a consequence: 

log y a = log|/a + log(l + (h/ a ), 
log(l + y a ) = log(l + y° a ) + log(l + y a Sy a ), y a = y^~ 0 - 

The source term is also expanded around its asymptotic counterpart: f a = f° + Sf a - 
Then equations (13.121) can be reformulated in terms of the 5y a functions as follows: 

log(l + Sy a ) = Sf a + Gab ~k log(l + y b Sy h ). (3.14) 

To dehne the iterative method, (13.141) are reformulated so that only 0(5yf) terms 
remain on the right hand side of the equations. Thus the equations are rewritten in 
the form: 


Sy a - Gab * Wvb) - 8 fa = Gab * [log(l + T/, Sy b ) - 34 Sy b ] - [log(l + 5y a ) - 5y a ]. 

(3.15) 

It can be seen that on the left hand side of (I3.15|) all the quantities are 0(Sy a ), while 
on the right hand side all the quantities are O^dy'l). This separation allows us to 
dehne an iterative solution. If Sy a s are small then the RHS is a small correction 

15 For repeated indexes summation is understood. 
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with respect to the LHS, this is why in an iterative solution the RHS can be simply 
taken at the value of the previous iteration. 

Let by a"' 1 the value of Sy a after the nth iteration, then byi” +1 ' ) can be determined 
from 5y ^ by solving a set of linear integral equations: 


5y^-G ab ^(y b 5yi n+1) )-5f a 


Gab * 


log(l + 34by;; )-y b 5y b 


[log(l + {»<">) - {»<")] . 


(3.16) 

Thus at each step of this iterative method a set of linear integral equations must 
be solved. Using the discretization method of the previous subsection, the problem 
reduces to solving a set of linear algebraic equations, which is a straightforward task 
in numerical mathematics. The very first (Oth) iteration starts from the asymptotic 
solution 5y a = 0 and it corresponds to the solution of the linearized equations, which 
in our case gives the Liischer-formula for the energy. 

This (first) method in a certain range of the coupling constant defined a numeri¬ 
cally convergent iteration to solve the equations for the ground state of our D-brane 
anti-D-brane problem, but beyond a certain value of g the method failed to converge 
anymore. This is why we worked out a second method, which proved to be much 
more efficient than the first one. This efficiency is manifested in two facts. First it 
converges much faster than the previous iterative method, second it gives convergent 
solutions to our equations even when the solution cannot be accepted as physical 
ont0- 


This second method can be described simply in words. Instead of defining an 
iteration as above, we simply take the discretized version of (j3.14j) . We consider 
it as a set of nonlinear algebraic equations. As a first step we solve the linearized 
discrete equations (i.e. (13.161) with RHS = 0) and starting from the solution of the 
linearized equations we solve the discrete nonlinear system by Newton- methods. 


16 Beyond a certain value of the coupling constant the equations in the form presented in section 
2. are not the right ones anymore, they should be corrected by some new source terms and 
quantization conditions, but even for the ’’wrong” equations the second method shows numerical 
convergence, giving unacceptable result. 

17 In MATHEMATIC A language it can be implemented by FindRoot [..., Method—>” Newton”]. 
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4 Numerical results 


In this section we summarize our numerical results. We solved numerically the 
equations for several integer values of the length parameter L. In this section we 
concentrate on the states with L > 2. The L = 1 special case is discussed in the 
next section. For the explanation of the numerical data we will mostly use the L = 2 
case as an example, because the critical point of this state is the closest one to zero, 
so it is enough to work with relatively small values of the coupling constant. This 
is important from the numerical point of view, since by increasing g the numerical 
method becomes more and more time consuming. 

First the parameters of the numerical method is discussed. There are three 
parameters in the nonlinear integral equations fl2.4[) - fl2.35D . The most important 
one is the coupling constant g , then there are two other parameters which allow us 
to formulate the equations according to our purposes. The two parameters are p 0 
and C, where po is a kind of ’’truncation index”, which tells us the node number 
starting from which the upper TBA equations are replaced by 577(4) NLIE variables 
(see figure m). The parameter C is a free parameter in the asymptotic solution 
for the upper SU(4) NLIE variables flC.6HC.20]) and it enters the equations such 
that the asymptotic solution around which the equations are formulated contain 
this parameter. From this discussion it is obvious that g is a physical parameter 
which means that the energy depends on it, while the other two parameters po and 
C correspond to different formulations of the same mathematical problem, so the 
energy does not depend on them. Thus the choice of these parameters is in our hand 
and we tried to choose such values for them which allows us numerical convergence 
in the widest range in g. For example the C = 0 choice is the best for numerical 
purposes since due to (12.221) a u —* — u symmetry arises in the S77(4) HNLIE 
variables minimizing the number of unknowns in the problem. Tuning p 0 might have 
two advantages. First, numerical experience shows that for large p 0 the Chebyshev 
coefficients of the unknowns entering the formula (12.101) for fl, decay faster, which 
allows for higher numerical precision. Second also from numerics we learn that with 
Po fixed at certain values of g non physical results are obtained from the numerical 
solution of the problem. This is a consequence of new local singularities entering the 
problem, but we still did not take them into account in the equations. These new 
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Figure 1: The pictorial representation of the E-system and the HNLIE structure with the choice 
Po = 4. 

singularities show up mostly in the 577(4) NLIE variables, thus by increasing the 
value of p 0 the appearance of such singularities in the equations can be postponed 
to higher values of g. 

We solved numerically our equations for different values of L and with various 
values of po and CL and in case the numerical result was physically acceptable for all 
Po and C we tried, it was also independent of these parameters within the numerical 
errors of the method. 

So far we discussed the parameters of the continuous integral equations and their 
role in the numerical solution. Now we turn to discuss the numerical parameters of 
the equations. The numerical parameters are artifacts of the numerical method, and 
they arise mostly from the discretization method described in section 3. We note that 
there is no cutoff parameter in our numerical method, neither in the integration range 
nor in the index of Y -functions. Everything is treated in an exact manner, the only 
source of numerical errors is the discretization of the unknowns and the convolutions. 
Here we give the most used subiuterval vectors of our numerical computations. 
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On each subinterval we used an 1%. = 10 order Chebyshev approximation. The 
subinterval vector A± of [0, 2g\ is given by the empirical formula: 


A + = 


A±,< if 9 < 2, 
A±,> if g > 2, 


(4.1) 


where the vectors in components take the form: 

2 gk 


A {k) = 

±’< [2g] + 1 ’ 


k — 1,[2g] + 1, 


A k ) _ 
A ±,> ~ 


1 3 11 

2 5 ,_ 4 ’ 2 5 ,_ 2' 2g ~ Z ,2g 


with v having vector components: 

2 g- 


Vj = l+ J 




3 = 1 ,-, 


29 - 1 
y 2 


(4.2) 

(4.3) 

(4.4) 


Here [...] stands for integer part. The set of subinterval vectors Aq of [0,H(<5)] 
could also be given by an appropriate empirical formula, but it would take such a 
complicated form, that it is better to write down the requirements from which it 
can be constructed. The requirements can be formulated in the language of the 
variable t G [0,-B(Q)]. The elements of the vector Aq divide the interval [0,H(Q)] 
into subintervals. Our requirements constrain the allowed length of the subintervals 
with respect their location within the whole interval [0, B(Q)]. The requirements 
are as follows: 


• The first clement of Aq is 

• The length of subintervals At in the range | < t < 2 is approximately 
At < |. 

• The length of subintervals in the range 2 < t < 3 is approximately |: At ~ 

• The length of subintervals in the range 3 < t < B(Q ) is approximately 1: 
At < 1. 

18 These requirements are based on numerical experiences with the choice h > 10. 
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In practice the length of the subintervals are slightly ’’squeezed” with respect to the 
conditions above to fill the full [0,£>(Q)] proper! \l 1 ’1. 

Finally, we note that for checking the numerical precision, we also did numeri¬ 
cal computations with 4 = 12,14,16 keeping the subintervals fixed and also with 
keeping 4 = 10, but doubling the number of subinterval points. 

Before turning to present the numerical results we would like to say a few words 
about the possible tests of the numerical results. Namely, how one can recognize a 
wrong result. This is also a very important point of the numerical method, since 
there are a lot of equations with very complicated kernels and it is easy to make 
mistakes during writing the code of the numerical solution. There are three basic 
things that we can check from the numerical results. 

The first check is dictated by the energy equation (12.121) . It is known that the 
energy starts at the first wrapping order (i.e. e~ L ) and this first order correction is 
exactly given by the Luscher formula [4]: 

OO °° 

= («) 
Q =1 o 

with Yq{u ) given explicitly in (ID. 21) . This quantity can be computed numerically 
with any digits of precision, so its value is known exactly at any values of g and 
L. The Liischer-formula (14. 5 1) corresponds to the linearized version of our equations 
(I2.4I) - (I2.35I) . this is why solving the linearized set of equations (which is the first 
step for the iterative solution) we should reproduce the numerical evaluation of 
(14.51) . This is a nontrivial check on the kernels, on the discretization method and on 
the equations themselves as well. In addition since this test is quantitative it can 
tell some information also on the numerical precision of the methocjfT 

This test can signal problems on solving the linearized problem. The remaining 
two tests can signal some discrepancies during the solution of the nonlinear problem. 

The second testing condition is that from the numerical solution Yp 0 _ 2 \vw must be 
real. This sound trivial, but it is not trivial at all. If one takes a look at the equation 
19 Not to have very small subintervals: At ^ 0.1. 

20 If one experiences that the numerical solution of the linearized problem agrees with the nu¬ 
merical value of (Id.51) within certain digits of precision, than the deviation from the Liischer result 
can be a good starting estimate to the numerical error. One cannot expect better accuracy, but 
the precision will not become much worse either. 
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(I2.32p of Y po _ 2 \ vw , one can recognize that there are complex quantities on the RHS 
which do not form conjugate pairs. So, the reality of the LHS is not guaranteed by 
the form of the equations, but it is guaranteed by the form of the solution. Thus 
the second testing condition is expressed by the inequality: 

|Imlog Y Po - 2 \vw\ A Numerical error, 10 -6 Numerical error 1CT 9 . (4.6) 

Here we wrote the typical numerical errors we had during the computations. 

The third test is based on the approximation scheme we use. One must check 
whether the Chebyshev coefficients of the unknowns decay as it is expected. From 
such a check the numerical precision of the method can be read off and it can shed 
light on some anomalous divergent behavior of the numerical solution. Thus it can 
indicate possible errors in the elimination of the divergent In u terms in 02.21) and 

The numerical results for the L — 2 case can be seen in figure [2] and table |T| 
In the table we show not only the energy E BTB a at different values of the coupling 
g, but the constant Sc, as well. The other columns of the table are related to the 
solution of the linearized equations; E BBBA and Sc^ are the energy and the global 
constant from the numerical solution of the linearized equations. A E BBBA stands 
for the deviation of E BBBA from the exact Liischer result. This quantity gives some 
information on the numerical accuracy of the method. Finally the column ’’number 
of nodes” tells us the cutoff index of the Liischer formula, which is necessary to 
get the Liischer energy with the precision given by AE bbba . This number is not 
equal to p 0 in our equations. For the L — 2 state, in case of 0 < g < 1.9 we used 
p 0 = 4, for 1.9 < g < 2.1 we used p 0 = 8, and in the range 2.1 < g < 2.14 we took 
Po = 12. Finally at g = 2.16 we used po = 26 to get acceptable numerical results. 
Then beyond this point we could not save our equations from the entrance of new 
singularities by increasing the value of po with a reasonable 0(10) amount. Because 
of this reason we could not get really close to the supposed critical point. There 
E B tba ~ —1; but we could reach only E BTB a ~ —0.7 at g = 2.16. Apart from this 
very embarrassing fact, some important features can be read off from the numerical 
data. First of all it can be seen that in the range g < 2.16 the energy is very slowly 
varying function of g, so there is no sign of any divergent behavior. What is more 
interesting is the behavior of the global constant Sc. It is negative and it decreases 
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faster and faster as g is increased. From the definition of 5c (12. 2 p it follows that all 
lQ-fnnctions are proportional to its exponent: Yq ~ £ = e Sc . The fast decrease of 
5c indicates that though Yq has worse and worse large u asymptotic by the increase 
of g. its global magnitude is actually decreasing. This remark can be understood 
from the TBA formulation of the energy. 


J BTBA 


oo 

£ 

0 = 1 ' 


tJtt' ogd + W). 


(4.7) 


Close to the critical point Ebtba is supposed to be finite [0] Ebtba ~ 1 — L, but 
naively the sum in the RHS of (I4.7|) would diverge due to the large Q terms. Since 
Yq is small for large Q, in leading ordecj the log(l + Yq) —> Yq replacement can be 
done: 


Q o 


JBTBA 


0=1 ' 


Zl£%W + Y Ql u))-( v l^y Q (u) +..., (4.8) 


Q—Q o q 


du dpQ 
27r du 


— 

Finite 


Diverges close to the critical point 


where Q 0 is an arbitrary index cutoff scale and Yq = £ Yq replacement was applied. 
Since £ is Q-independent all the dangerous Q dependence is still in Yq. In (14. 8 p 
approaching to the critical point the second sum starts to diverge, and the global 
multiplicative factor £ must tend to zero in order to ensure the finiteness of both 
sides of the equation. Our numerical data seems to support this picture. Namely 
5c -A- —oo as going closer and closer to the critical point. 

In [5] from Y -system arguments the large Q behavior of Yq was also estimated 
by the formula: 

y e(“> - «») 7- Asw y e(“). « = ^ ( 4 -9) 

(“ 2 + f ) 

where 5c is defined after (12.2p in section 2. (14.9p is a very important formula, 

because it plays crucial role in the analytical determination of the critical point. 
Since the numerical solution of the HNLIE equations of section 2. does not require 
the introduction of any index cutoff, it takes into account the contributions of all 

21 For large Q. 
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Ebtba 




Figure 2: Ebtba (on the left) and 5 c (on the right) as functions of g for the L = 2 state. 


9 

Ebtba 

5c 

r(°) 

^ BTBA 

Sc® 

Ar(°) 

BTBA 

number of nodes 

1.6 

-0.175553 

-0.844383 

-0.185898 

-0.893355 

3.7 

• 10~ 6 

25 

1.7 

-0.2271599 

-1.17751 

-0.24183601 

-1.21357 

4.9 

• 10~ 6 

25 

1.75 

-0.25693719 

-1.36622 

-0.2738668 

-1.40256 

5.7 

• 10~ 6 

27 

1.80 

-0.2897776 

-1.58077 

-0.3088130 

-1.61278 

6.3 

• 10~ 6 

28 

1.90 

-0.366494169 

-2.10766 

-0.38810198 

-2.10321 

7.9 

• 10" 6 

30 

1.92 

-0.38393979 

-2.23237 

-0.40555472 

-2.21339 

8.2 

• 10~ 6 

30 

1.94 

-0.402255118 

-2.36573 

-0.4235649 

-2.32785 

8.7 

• 10“ 6 

30 

1.96 

-0.42147149 

-2.50781 

-0.4421440 

-2.44671 

9.0 

• 10~ 6 

30 

2.00 

-0.46303978 

-2.82377 

-0.4810544 

-2.69809 

9.9 

• 10" 6 

31 

2.02 

-0.48564199 

-3.00085 

-0.5014098 

-2.83086 

9.9 

• 10" 6 

32 

2.04 

-0.50966430 

-3.19333 

-0.52237993 

-2.96847 

1.0 

• 10~ 5 

32 

2.06 

-0.53532776 

-3.40422 

-0.5439774 

-3.11107 

1.0 

• 10“ 5 

32 

2.08 

-0.56291307 

-3.63744 

-0.566214 

-3.25878 

1.0 

• 10~ 5 

33 

2.10 

-0.592805 

-3.89861 

-0.589106 

-3.41179 

8.9 

• 10~ 6 

35 

2.12 

-0.625515 

-4.19506 

-0.612655 

-3.56999 

1.2 

• 10~ 5 

34 

2.14 

-0.661868 

-4.54055 

-0.636888 

-3.73339 

9.6 

• 10~ 6 

36 

2.16 

-0.7031687 

-4.956683 

-0.661809 

-3.90338 

7.0 

• 10~ 6 

39 


Table 1: Numerical data for the L = 2 state. 
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the Y -functions of the infinite Y -system. This makes it possible to test numerically 
the correctness of the large Q estimate (14.91) . In case (14.91) holds, it implies that 
8 In jjq = In ijq — In T/q tends to zero as 1/Q for large Q. In figure El the numerical 
demonstration of this statement can be seen. The plotted functions are defined by 
the formula: 


SFcj{t) = 


5\ny Q (x q (B(Q) - t)) if t > 0, 
dlny Q (-x Q (B(Q)+t)) if t < 0, 


(4.10) 


where x Q (t) = Q 


Lj , B(Q ) = ^ + 1. The plots of figured! are based on 
the numerical computation with p 0 = 26 at g = 2.16. Figured) nicely demonstrates 
the expected 1/Q behavior of the functions 5Tq. 



Figure 3: Large Q behavior of SYq from numerical data at g = 2.16 with po = 26. 

For the L = 2 state beyond g = 2.16 the numerical solution of the discretized 
problem did not give physically acceptable results. To get some insight into the 
source of the problems, at g = 2.18 we plotted the imaginary part of the LHS of 
the last equation in (12.18j) . namely Imlog(l + SbQ at u — x Po (B(p 0 ) — t). Figure 
[H shows that there is a jump of 27 t, when t is close to B(p 0 ). (I.e. large u. j^l 

22 Here the sampling points are connected according to the Clrebyshev approximation. This is 
why the jump of the logarithm is not ’’sharp”. 
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This fact shows us that the equations we solved numerically are not the right ones 
anymore. Something is missing from the equations. Either a special object gang 
or some other local singularities of the T- and Q -functions of the problem, which 
enter those strips of the complex plane, which are relevant in the derivation of the 
ffNLIE equations. 

The numerical data for the L — 3 and L — 4 states are given by table El and El 
Also in case of these states the appearance of new singularities obstacled us to get 
close to the critical point in the framework of the HNLIE technique. 



Figure 4: The anomalous behavior of Imlog(l + Sbe) at g = 2.18 and po = 26. 


5 Comments on the L = 1 case 

The L = 1 ground state is mysterious, since so far the anomalous dimension of 
this state could not be determined even for small g either from field theory or from 
integrability considerations [4] . ffere we concentrate on the integrability side. There 
the boundary Liischer formula HUES] diverges for this state [4]. For generic L the 
Liischer formula is simply the expansion of the TBA energy formula around the 
asymptotic solution with the replacement: log(l + Yq) —> Yq. For small coupling it 
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9 

Ebtba 

5c 

^ BTBA 


AP(°) 

L ^ Cj BTBA 

2.2 

-0.114591 

-0.62945 

-0.12907 

-0.711869 

8.4 

h- 1 

o 

1 

-4 

2.6 

-0.21909 

-1.33443 

-0.267823 

-1.641000 

3.0 

o 

1 

-4 

2.8 

-0.286833 

-1.82583 

-0.366547 

-2.34855 

2.1 

• 10~ 6 

3.0 

-0.365866 

-2.42503 

-0.488968 

-3.26271 

2.9 

• 10“ 6 

3.2 

-0.457294 

-3.14677 

-0.638504 

-4.42131 

8.7 

• 10“ 6 

3.4 

-0.56282 

-4.01232 

-0.818842 

-5.86636 

1.1 

• 10" 5 

3.6 

-0.685108 

-5.05271 

-1.03391 

-7.64329 

1.5 

• 10~ 5 


Table 2: Numerical data for the L — 3 state. 
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Ebtba 

5c 

^ BTBA 

5c^ 

ApM 

Lxrj BTBA 

2.6 

-0.0716174 

-0.427755 

-0.0793412 

-0.476413 

6.0 

• 10” 7 

2.8 

-0.0975242 

-0.607523 

-0.111564 

-0.699788 

1.1 

• 10“ 8 

3.0 

-0.128116 

-0.829046 

-0.151888 

-0.991352 

4.2 

• 10“ 8 

3.2 

-0.163439 

-1.0949 

-0.201422 

-1.36341 

9.2 

• 10” 8 

3.4 

-0.203514 

-1.40733 

-0.261362 

-1.82949 

9.8 

o 

1 

00 

3.6 

-0.24835 

-1.76832 

-0.332987 

-2.40437 

7.0 

•10" 7 

4.0 

-0.35239 

-2.64384 

-0.51686 

-3.94631 

8.6 

o 

1 

-4 

4.2 

-0.411691 

-3.16247 

-0.632106 

-4.9496 

2.7 

• 10~ 6 

4.4 

-0.545354 

-4.3733 

-0.917326 

-7.52081 

1.3 

• 10~ 6 


Table 3: Numerical data for the L = 4 state. 
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takes the form [4J: 


A E(L) 



du dpQ 
27r du 





I) 


4 L 


4L — 1 


4L 

2L 


C(4L — 3) + 0(g 2 ) 


(5.1) 

This small coupling expression diverges for L = 1, since this point sits exactly 
on the pole of the ( -function. As for the origin of this divergence; in (15. lj) the 
individual integrals are convergent, but their sum for Q causes the divergence. In 
|4j it was argued that also for any larger L the TBA energy formula would diverge 
beyond a certain critical value of the coupling: g c {L). Assuming that the energy is 
a monotonously decreasing function of g, which is supported by numerical results, 
this critical point can be expressed clearly in terms of the energy by the criterion: 


E C (L) = E(g c (L)) = 1 - L. 


(5.2) 


In [4] this point was interpreted as a turning point where the energy becomes imag¬ 
inary and as a physical consequence the ground state becomes tachyonic. For the 
L — 1 state the critical point is right at g = 0 assuming that for small g the energy 
is also small. 

Now let us turn our attention to the HNLIE description of the problem detailed 
in section 2. Here there are no infinite sums and even for L = 1 all the convolutions 


of the integral equations seem to converger 3 !. For the first sight there is no sign of 
any problem in the HNLIE description and it seems that only the TBA description 
is inappropriate to treat the L — 1 case. But unfortunately this is not the case. 

We can write down the discretized integral equations for the L = 1 case as well, 
and using the Newton-method, we can solve them for small values of the couplinaF 3 !. 
We always get some numerical solution for the discretized problem, but it turns out 
that the Chebyshev coefficients of the unknowns, which correspond to the large u 
subinterval do not form a decaying series. This phenomenon is a typical sign of some 
weak (probably logarithmic) large u divergence of the unknowns. If one increases 
the number of subintervals and sampling points the situation remains the same. The 


23 If we assume that large u behavior of the unknown functions, which was used to derive the 
BTBA equations from discontinuity relations and Y-system. 

24 Typically g ~ 10 _1 . 
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conclusion is that we can solve the discretized problem, but the solution cannot be 
interpreted as the discretely approximated version of the continuous solution of our 
integral equations. In other words the continuous HNLIE equations have no solution 
for L — 1. 

In order to get some analytical insight why the solutions become diverging at 
large u let us consider the TBA formulation of the problem (p 0 —> oo in HNLIE). It 
is known [4] that the TBA energy comes from the coefficient of the most divergent 
log \u\ term in the large u expansion of log Yq: 

log Y q (u) = —4 (L + Ebtba ) log \u\ + 0(1). (5.3) 


The Ebtba term originates from the RHS of the TBA equations for log Yq from 

OO 

the convolution term log(l + Yq>) * by exploiting the large u expansion of 

Q '=l 

the kernel: K®^(v,u) = log \u\ + 0(1). has better large Q' behavior 

than that of -^n-, since it behaves like 1 /Q'. As a consequence contrary to the 
energy formula, the sum of dressing convolutions is convergent indeed. Thus one 
might think that for L — 1 the problem emerges, because for the derivation of 
the energy formula we expanded the sum of dressing convolutions term by term 
for large u. This is why instead of this usual procedure, we consider the sum of 
dressing convolutions itself, compute it and then at the end of the computation we 
take the large u expansion. This procedure is carried out in the small coupling limit. 
We need the leading order small coupling expression of the dressing kernel in the 
mirror-mirror channel^!: 




1 

2n 


lM 1 + y-q« I )+^( 1 + y + q« I 


/ 11 i Q' + Q i • w a 
-V ( 1 + —2 — + l 2 _ Ul ' 


<P I 1 + _ jl(« 2 - «,) 

(5.4) 


Then the formula, the large u expansion of which accounts for the small coupling 
expanded energy, is given by: 



Ch y 


47T 


o,(L= 1) 
Q' 


v)K 


Q'Q,( o), 

«i(2) 1 


(5.5) 


+ • 


25 Here we use the rapidity convention where the branch points are at ±2 g. 








where Yq^^^u) = g 4 ^ 2 +q 2 ) 3 is the leading small coupling expression of (ID. 211 at 
L = 1. The second derivative of 0(u,Q) can be computed explicitly by simple 
Fourier space technique. We take the Fourier form of each functions under integra¬ 
tion, the convolution is the product of the individual Fourier transforms, the sum 
for Q' can be easily done in Fourier space and at the end of the process everything 
is transformed back to the u space. In such a manner one gets a bulky, but explicit 
expression for J^(9( u , Q ), which we do not present here, only its large u expansion: 

Y^O(u,Q) = ig 4 — 7 ^ + 21112 ^) \ - 8g 4 + O(^). (5.6) 

au- \2 J u z u z ir 

Integrating twice the large u expansion at small coupling becomes: 

0(u 1 Q)=Ag 4 Q + 27 s-In4^ logu + 4/(logu) 2 + ... (5.7) 

From (15.71) it is obvious why the naive Liischer energy formula diverged. Because 
the leading order large u term is not the expected ~ log |w|, but ~ (logn) 2 . This is 
the key point of the problem, since in this case after this first iteration Yq acquires 
an unwanted type of large u term, which makes Yq divergent for large u: 

Y Q (U ) ~ u (-±L+ig*(\+2 lE -ln4)+...) e 4 fl 4 (l 0gU ) 2 +..._ ( 5 . 8 ) 

This large u divergence contradicts to what was assumed about the large u behavior 
of Yq at the derivation of the integral equations, since it was supposed to decay. In 
this example we have shown in the small coupling limit, that during the iterative 
solution of the BTBA equations, log Yq acquires an extra ~ (log|w|) 2 behavior at 
infinity, which made Yq an exploding function at infinity. This means that the 
iterative solution of the TBA equations leaves the class of physically acceptable 
solutions. 

One might ask the question, whether it is possible to keep somehow the qualita¬ 
tive large u behaviors that we assumed at the derivation of the equations? Here we 
sketch a possible idea for small coupling to the L — 1 case. 

Let us assume that we managed to modify the TBA equations, such that all 
Y -functions have the large u behavior we want. Since most of the TBA equations 
reflect the structure of the Y -system functional equations we expect to modify only 
those equations which are affected by also the discontinuity relations. It follows, 
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that for large Q, the formula for the estimate for Yq (j4.9j) remains the same. Now, 
we assume that for small g the energy is also small and take the simultaneous small 
g and small energy expansion of the RHS of the TBA energy formula (14.711 . In 
leading order the large Q terms will dominate: 


Ej 

Q= i_“ 


du - f u 

4^ Ms 


~'E J 

Q= l_ r 


du / 2\ 2L inn 2 _ — _ 

^ \9 ) h: ,2 _|_ Q2^2(L+E B tba)+^ 




4.9 


4L 


4(L + E BTBA )\ 


E 


^(L+Ebtba) 4(L + E btba ) - 1 \2(L + E btba ) J ^ Q( 4 R+Estba)-3) 
{ ' 4941 ^ 4(i + + £ btba ) - 3) 


24(L+s srsA ) 4(L + E BTBj 4 ) — 1 \2(L + E btba ) 
I.9 4 




BTBA 


+ 0(M 4 ), 


(5.9) 


where Eq denotes the large Q estimate (j4.9j) of Yq, £ = ^ g' iEBTBA as a consequence 
of the u —y u/g change of variables and the pole term in E BB ba comes from the 
pole of the ((-function. I 11 om 1 HNLIE approach the energy E BB ba and the constant 
Sc are parts of the equations which means that they are not simply expressed by 
explicit formulas based on the solution of the equations, but must me obtained by 
solving the set of non-trivially entangled equations. In this sense (15 .9j) defines an 
equation for E BB ba for small g. Its leading order solution is: 


Ebtba — g 2 



(5.10) 


If £ > 0 then E BB ba becomes imaginary as it would be expected from string-theory 
expectations [4j. To decide the sign of £, the equation (12. 131) has to be analyzed in 
the context of the small g and E BB ba expansion. It turns out that £ is positive and 
0(1) for small g, so according to (15. 10|) E B tba is imaginary. Another remarkable 
fact is that according to (15. 101) E B tba starts at 0(g 2 ) instead of the 0(g 4 ) prediction 
of the boundary Liischer formula (15.11b This might be another explanation why the 
coefficient of g 4 diverges in the Liischer formula for the L = 1 case. Finally, we note 
that in the small g and E BB ba expansion of the L = 1 state, the energy is pure 
imaginary only at leading order in g , but in higher orders it acquires real part as 
well. 
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For the first sight, it might seem that without modifying the equations one 
immediately gets imaginary energy when going through the critical point. But, the 
situation is a bit more subtle. There is a hidden tacit modification of the equations. 
This is realized in (15. 9ft by the replacement: 

OO ^ 

^2 Q4(L+E BT ba)-3 C( 4 ( R + E BTBA ) - 3). 

For the L — 1 case it is an identity for R q{Ebtba) > 0, but for R s(Ebtba) < 0 it 
is not an identity anymore, but a nontrivial analytical continuation in Ebtba- 

Such an analytical continuation would require the exact determination of com¬ 
plicated sums of convolutions of the TBA equations as functions of the energy. Since 
this does not seem to be feasible in practice, we give such an alternative modification 
of the TBA equations which preserves the infinite sum structure of the equations, 
but the sums will converge everywhere for R q^Ebtba) > —L except at the critical 
value Ear = 1 — L. 

The basic idea of the modification comes from the sum representations of the 
^-function. The usual one converges for Re(s) > 1: 

oo 1 

co) = Re 0) > o (5.ii) 

Q =l ^ 

but there is another representation which converges for Re(s) > 0: 

<w = rh£(«TiF-^r)- Re(s)>a (5 ' 12) 

Then the original TBA equations are modified through their infinite sums by the 
replacements: 

OO ^ OO 

^ Lq*Eq -» - - {Q ■ {Lq +1 -k JCq+i) — (Q — s E ) ■ ( L q -k ICq)} , (5.13) 

Q=i _ 1 0=1 

where s E = 4(L + Ebtba ) — 3. Taking into account the large Q behavior of all 
Yq functions and all the kernels of the inhnite sums of the TBA equations, the 
new representation will converge for R g(Ebtba) > ~L. This slight modification of 
the TBA equations might make it possible to go beyond the critical point and get 
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solution of the TBA equations with large u asymptotics being in accordance with 
the ones used for the derivation of the equations. 

The conclusion of this heuristic argument is that to keep the expected^! qualita¬ 
tive large u behavior a nontrivial modification of the TBA equations must be carried 
out, which might lead to complex energies. 

6 Summary and conclusions 

In this paper we studied the ground state energy of a pair of open strings stretching 
between a coincident D3-brane anti-.D3-brane pair in S 5 of AdS§ x S 5 . The main 
motivation for the study is that string-theory predicts that the ground state of such 
a configuration becomes tachyonic for large values of the ’t Hooft coupling [4j. 

In [4] it was shown that the usual integrability based BTBA approach always 
give real energies for the ground state and it breaks down at latest when the energy 
gets close to the critical value: E C (L) = 1 — L. This point was interpreted in jl] as 
a transition point where the ground state becomes tachyonic. 

Approaching this critical point the contribution of all the Y -functions of the 
BTBA becomes quantitatively relevant, thus the numerical solution of the trun¬ 
cated BTBA equations cannot give accurate results close to the critical point. To 
resolve this difficulty and get more accurate numerical results we transformed the 
previously proposed BTBA equations into finite component HNLIE equations. The 
HNLIE equations were solved at different values of g and L and the numerical results 
confirmed the earlier BTBA data. 

During the numerical solution of the HNLIE equations the usual iterative meth¬ 
ods failed to converge, this is why we worked out two numerical methods to reach 
convergence. The most effective one is, if one transforms the integral equations into 
discrete nonlinear algebraic equations and solves them by Newton-method. The 
power of this method is demonstrated by the fact that it gives convergent results 
even if the numerical solution is not physically acceptable. 

Unfortunately, in our numerical studies we could not get very close to the critical 
point, because new singularities entered the HNLIE equations taking into account of 

26 This primarily means that log Yq ~ log |u| for large u, while other T-functions tend to constant. 
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which would have required an enormous amount of additional work. Nevertheless, 
in the range where we could get physically acceptable results, the precision of the 
HNLIE data were higher than those of BTBA and the HNLIE approach could give 
a deeper understanding of the problem. 

For the ground state of the L — 1 state the critical point is right at g = 0 and 
neither perturbative field theory computations nor the boundary Liischer formula 
could provide a finite quantitative answer to the anomalous dimension. Even in this 
special case the numerical solution of the HNLIE equations was possible. The results 
showed that without an appropriate modification of the equations, they cannot give 
physically acceptable results. In this case, it means that the solution of the dicretized 
problem cannot be considered as a discretized solution of the continuous nonlinear 
integral equations. Moreover the large rapidity behavior of the numerical solution 
is incompatible with the one assumed for the derivation of the equations. This 
phenomenon is analytically analyzed in the framework of BTBA and an idea is 
sketched to preserve the expected large rapidity behavior of the unknowns. This 
method is based on an appropriate modification of the TBA equations which would 
lead to complex energies beyond the critical point. 

Hopefully the L — 1 case at g — 0 could be treated analytically in the framework 
of the quantum spectral curve method mmm, solving the mystery of this state 
in the context of integrability. 
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A Notations, kinematical variables, kernels 


Throughout the paper we use the basic notations and TBA kernels of ref. |4T], 
which we summarize below. For any function /, we denote f ± {u ) = f(u ± and in 
general f^ ±a ‘(u) = f(u±^a), where the relation between g and the’t Hooft coupling 
A is given by A = 47r 2 g 2 . Most of the kernels and also the asymptotic solutions of 
the HNLIE-system are expressed in terms of the function x(u): 


x(u ) = —{u — iV4 — u 2 ), Im x(u) < 0, 


(A.l) 


which maps the u-plane with cuts [—oo, —2] U [2, oo] onto the physical region of the 
mirror theory, and in terms of the function x s (u) 


x s (u) 




x s (u)| > 1 , 


(A.2) 


which maps the w-plane with the cut [—2, 2] onto the physical region of the string 
theory. Both functions satisfy the identity x(u) + yyy = u and they are related by 
the x(u) = x s (u), and x(u) = 1 /x s (u) relations on the lower and upper half planes 
of the complex plane respectively. 

The momentum p ® and the energy Sq of a mirror Q-particle are expressed in 
terms of x(u) as follows: 

i / i 

u - Q) — gxlu H— 

9 9 

Two different types of convolutions appear in the HNLIE equations. These are: 

/ oo /*2 

du f(u) IC(u,v) , / */C(u) = / du f(u) K,{u,v). 

-oo J —2 

The kernels and kernel vectors entering the HNLIE equations can be grouped into 
two sets. The kernels from the first group are functions of only the difference of 
the rapidities, thus actually they depend on a single variable. The other group of 
kernels composed of those, which are not of difference type. 


p Q (u) = gx( 


X \ 11 - ~~\cJ ) 

Q) + tQ, £ Q (u) = \oe x(u+ l Qy (A.3) 
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We start with listing the kernels depending on a single variable: 


s( “ ) = 2^ logT " ( “ ) = 


T(u) = tanh[ f“! ■ 


2 

gQ 


Kq ^ 2m du l ° gSQ ^ 7 r Q 2 + 3 2 « 2 

1 d 


, $,(«) = 


U 


iQ 


u+iQ 


9 

M—l 


K mn (u) — -y— — log S M n(u) — K m+n (u) + K n _ m {u ) + 2 E KN-M+2j{u) , 
m u , =i 

M -1 

Smn{u) = Sm+n{u)Sn-m{u) n SN-M+2j(u) 2 — Snm(u ) . (A.4) 


J=1 


The fundamental building block of kernels which are not of difference type is: 


1 d 


1 V^T 


K(u, v) = — — log S(u, v) = — - 

2m du 2m 1/4 — u 2 u — v 


, S{u,v) = 


x(u) — x(v) 


. (A.5) 


x(u)x(v) — 1 

Using the kernels K(u,v) and Kq(u — v ) it is possible to define a series of kernels 
which are connected to the fermionic Y±- functions. They are: 


X X 

K q Ju,v) = K(u - Q,v) - K(u + -Q,v), 

9 9 


KS y (u, v) = 4 


u-v) ± K Qy (u,v)j 


and 


X X 

KyQ (■ u , v) = K(u,v + —Q) — K(u, v - -Q ), 
Kl Q (u, v ) = i (KyQ(u, V ) =F Kq(u - vfj . 


(A.6) 

(A.7) 

(A.8) 

(A.9) 


Further important kernels entering the Y± related TBA-type equations are defined 
as follows: 


K° m (u,v) = 


1 d 
2iri du 


log Sg M (u,v), 


x(u — i—) — x(v + i— ) x(u — i—) — x(v — i—) x(u + i- 

Sg M (u,v) = { 9> ^ 9> ^ 9> 1 9> 1 9 


x(u + l 


Q' 


x(v + *y) x(u + fy 


x(v — i—) x(u — i®-' 

v 9 ' v 9 ' 


11 — v — >- 4 / — 2 /) 


(A.10) 
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The kernels entering the right hand sides of the equation (12.91) for Y\ are 



x 


x(u — i—) — x(v + i—) x(u + i—) — x(v + i—) 
x(u — ij) — x(v — iy) x(u + ij) — x(v — i y) 

Y^u- V - ~ g (M -Q + 2j) 

U-V + i(Af -Q + 2j) ’ 


x(v — i—) 

+ f—) 

v 9 y 

(A.ll) 


and the dressing-phase related kernel K®^(u,v), which is built from the s((2) S- 
matrix of the model [49]. ft is of the form 

s S^)( u , v ) = S QM {u - n) _1 Y qm (u, v)~ 2 , (A.12) 

where is the improved dressing factor [50] • The corresponding st(2) and dress¬ 
ing kernels are defined in the usual way 

K Sp)( u > v) = iog S^{u, v), Kq M (u, = 1o S ^qm(u, v ). (A.13) 

Explicit expressions for the improved dressing factors Y,qm(u,v) can be found in 
section 6 of ref. [50]. Here for our numerical computations we used the single 
integral representation given in m 

Finally we mention that along the lines of [42] in the derivation of the formula 
(12.1 Op for f1 (/Cq), it was exploited that all the necessary kernels: 

K q , K Qy , A«\s * A-O-w, K„ U K% 2) satisfy the identity: 

IC Q - s -k )Cq-i - s */Cq+i = 5K,q = 0, for Q > 3. (A.14) 


B Kernel matrices of the vertical HNLIE part 


In this appendix the kernel matrices appearing in the upper HNLIE part of our 
equations (12.18112. 19p are presented. Here the kernel matrices are different compared 
to those published in [38]. The difference comes simply from a reformulation the 
equations in the language of new unknown functions. In |3B] the unknowns are 6 
6-type functions: 


b old _ r fc (3)[ 7 i] 6 (3)foa] 6 (3)[7 s] & (2)[-l+ 7 4] & (2)[-1+7 s] & (l)[-2+7e]| 


(B.l) 
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and 6 d-type functions: 


jold _ r j(3)[»?i] t(3)[)? 2 ] j(3)[j? 3 ] t(2)[jj 4 ] j(2)[jj 5 ] jp/faeh 

it — 1 U 1,S ) a 2,s 5 u 3,s > u l,s ? u 2,s ; U 1 ,s h 


(B.2) 


with shift vectors 7 and rj given by 


We recognized that the kernels 


become simpler if we formulate the equations in terms of the unknowns: 

i={^ [7 t>)/4!’ tal ,C ta ^/C hI+7 ‘ 1 .'>S hl+7 ‘ l .'> ( M h2+7 ’ 1 }. (B-3) 

and 




(B-4) 


where 7 / = ±1 is a global sign factor and s = po, if one adopts the notation of 
for the HNLIE equations ()2.18H2.19p . Another advantage of using the variables 
is that they are either 0(1) or exponentially small for large volumes. 



For the sake of simplicity, here we give the form of the kernels of (12.18H2.19p 
before the application of the contour shifts (j2.20H2.2ip . The kernels of the equations 
can be obtained from these by simply shifting their arguments according to the 
formulas below: 


GbB(u)ab — KbB Ul-(7a — 7b) 

9 / ab 


G bD {u) ab = K bD [u + - (j a - 7/6) 

9 / ab 


G ab I^cLB ( U T ('7/a 7 b) 

9 / ab 


G dD {u) ab = K dD (u + -(r] a - 7/ b ) 

9 y ab 


a,b = 1,..., 6 
a, b = 1,..., 6 
a, b = 1,..., 6 
a, b = 1,..., 6. 


(B.5) 


The kernel matrices can be expressed by the functions as follow: 


s 


N 9 igu, , igu. . ,1 igu . ,.1 igu . , . 

G(u) = <t/’( 1 H )+'0(l 4 ) — ^(2 ^ 4 ) — ^(2 ^ 4 ^ ’ (B'6) 


/ \ 9 1 

s(u) = - 


4 cosh(=f) 


(B.7) 

(B.8) 
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and they take the form: 


K h R 


( 


KbD = 


V 

/ 


KdB = 


V 


K(io = 



/ 

G 

0 

G-s+ 

—s 

s + -G 

0 



0 

l 

0 

l — s~ 

s 

0 



G-s~ 

0 

G 

0 

s~-G 

0 



s 

l - s + 

0 

l 

0 

s 



S — G 

—s 

s+ -G 

0 

G 

s~ 


V 

0 

0 

0 

—s 

s + 

G 


-G 

0 

s+-G 

G-s + 0 

0 

0 

-/ 

0 

—s s + — l 

0 

s~-G 

0 

-G 

G — s~ s 

0 

—s 

s + - l 

0 

0 s + -l 

s 

G - s~ 

s 

G-s+ 

1 

Q 

o 

s~ 

0 

0 

0 

s —s 

-G— 


-G 

0 

s + -G 

s 

G — s + 

0 

0 

-l 

0 

s~ — l 

—s 

0 

s~-G 

0 

-G 

0 

G-s- 

0 

G — s~ 

s 

G — s + 

0 

-G 

s + 

0 

s~ — l 

—s 

s - — l 

0 

s 

0 

0 

0 

—s 

s + 
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/ 

G 

0 

G-s+ 

s+ -G 

0 

0 
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0 

/ 

0 

s 

/ - s + 

0 



G-s- 

0 

G 

s~-G 

—s 

0 



s~ -G 

—s 

s+ -G 

G 

0 

s + 



0 

l — s~ 

s 

0 

/ 

s 


V 

0 

0 

0 


—s 

G 

/ 


(B.9) 


(B.10) 


(B.ll) 


(B.12) 
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C Asymptotic solutions of the vertical HNLIE 


In this section along the lines of [3Bj the asymptotic solutions of the upper 517(4) 
NLIE variables are presented . In the asymptotic limit the T-hook of AdS/CFT 
splits into two 5I7(2|2) fat-hooks. The basic building blocks of the asymptotic 
solution are the nine Q -functions corresponding to the left and right 517(212) fat- 
hooks. Due to the left-right symmetry of the Y -system it is enough to give the right 
Q-functions. They can be derived from the asymptotic solution of the Y-functions 
given in f4]. They take the form: 

Q( 1,2 \u) = g Aqu a{u) u + , & ’ \u) = 4 gu ++ , 

Q {1 ' 1] (u) = q u a{u ), Q m {u) = 

Q (1 ’ 0) (w) = Aq n a(u), Q m {u) = 1, 

(C.l) 


Q 


Q {2 


( 2 ’ 2 ) <u) = 


= 


= 922 , 

2 922 


g (2 ’ 0) («) = 


gAu 
4 q 22 u~ 
gu~ u~ 


where qn,q 22 and A are arbitrary constants which cancel from the final form of 
the asymptotic NLIE variables. Furthermore cr(u) = e~^~ to satisfy the recursion 
= — 1. The further building blocks of the asymptotic solution are as follows! 2 ^! 


77 , = 


4 (-l) s SM 




and 


A°{u) = 


4 u 


gu + u 


B°{u)=4gu, p\u) = 


2 cr {u) 

A : 


7°(u) = 


The solution of the recursions 


w°~ - w° + = 


A° 


o/O-t-^wO— 5 


y° + - y° = 


B° 


are as follows: 

, iA 2 e~ ngu 


w u = 


((gu) 2 + w c ), 


P°P°" ’ 

i A 2 e~ n9U 


(C.2) 


2 a{u) 

q Au 
(C.3) 


(C.4) 


y w = 


{{gu) 2 + w c - iC ), 
(C.5) 
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Here for correspondence we use the same letters for the names of different unknowns as in 


ESI- 
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where w c and C are arbitrary constants. Using the building blocks listed above, the 
asymptotic form of the upper £77(4) NLIE functions can be determined [38j and 
take the form: 


b[ 3 !°(u) 


b^°( u ) = — 

b 3,s W w [_ a] > 

p( 3 )°/..\ ( S + l ) U 

B3 :S ( U ) ~ r _i 


b?i°l 


u = 


m 


Bf]°(u) = -&- s \u) 


Ag su 


1 = 77" w = 
1 ,( 2 )°- 


s 1(m) 

Ag su + C 1 


C (“) = 


C + Agsu~ 
Agu^~ s } ’ 

C + 4 g (s — 1 ) u 
Agu^~ s ^ 


B, 


(2)o-, 


2 ,s 


n = 


C + 4 <7 (s — 1 ) u 


Agw- 


(C. 6 ) 

(C.7) 

(C. 8 ) 

(C.9) 

,(C.10) 


g + 4 ^7~ , (c. 11 ) 

4(yf m[ _S 1 






4 3 ]°(«) = 4 3 i » = 

= <» = (s+1) “ 
45 ° («) = b^ s o (u) = 


uH ’ 

1 ( gu W) 2 

01 s ! (u) Agsu 


(C.12) 

(C.13) 

(C.14) 


43» = 


C + Agsu + 


Agu M 

47 w = d?1°(u) = 

^’ s v ' 2 ’ s v ' 45-sn + C 

C + 4p(s-l)« (i )o 


=- c+4 / g ( ;,; 1)M . (c.i5) 


4!i°(«) = 


where s is the ’’cutoff index” where the TBA —> HNLIE replacements started 
furthermore for any index distribution B° and D° stand for 1 + 6 ° and 1 + d 
respectively. 


C + Ag su + 


(C.16) 

(C.17) 


</>(«) = 


x\u 


,2 L 


gu 


(C.18) 


and C is the arbitrary constant that does not cancel from the formula for the HNLIE 
variables. The asymptotic solution for the six b- and d-type NLIE-functions of the 


29 In section 2. it is denoted by po, here the notation s is kept to fit to formulas of [3Bj. 
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system can be obtained from the Backlnnd functions above by appropriately shifting 
their arguments: 


b° = m = >C 1 ' 

v(3)o[?7i] / t(3)o[? 72] 7(3)o[r/ 3 ] ,(2)o[r? 4 ] ,(2)o[r? 5 ] ,(l)o[»? 6 ] 


}, 

(C.19) 
(C.20) 

Finally we note that the C — 0 choice implies a 


d° = k} = J »»^7. J »>. 

with the shifts given in (12.20112. 2 
symmetry relation between the b- and d-type variables. Let 9ft the 6 by 6 matrix: 

/ 0 0 1 0 0 0 \ 

0 1 0 0 0 0 

1 0 0 0 0 0 

0 0 0 0 1 0 

0 0 0 1 0 0 

0 0 0 0 0 1 


m = 


V 


(C.21) 




Then at C = 0 the b° and d° vectors satisfy the relations as follows: 


d°0) = 9ftb°(—n), b°(—n) = b°*(u), (C.22) 

d°(n) = 9Jtb°*(w), d°(—n) = d°*( M ), (C.23) 

where * denotes complex conjugation. I 11 our numerical studies we mostly use the 
C = 0 asymptotic solution to setup the equations to solve. In this case the exact 
equations guarantee the fulfillment of (I2.22H . which reduces to 6 the number of 
independent complex functions of the upper NLIE part. 


D Asymptotic solutions of the F-system and the 
horizontal SU( 2)-type HNLIE 

This appendix is devoted to give the asymptotic solution for the K-functions and 
the variables of the horizontal SU( 2) NLIE. The asymptotic form of the Y -functions 
can be read off from the asymptotic T-functions in [4j. They take the form: 

2 ) g 2 u 2 


Y m\vw( U ) = 


m[m 


(■m + l) 2 + g 2 u 2 ' 


m = 1 , 2 , 


(D.l) 
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1 



(«) 


1 QQ 2 g 2 u 2 
g 2 u 2 + Q 2 ' 


<9 = 1,2,... 


(D.2) 


Y-(u) 

Yi\M 


Y?(u) = 


g 2 u 2 


2 + g 2 u 2 ' 


g 2 u 2 (19 + 3 g 2 u 2 ) 
(1 + g 2 u 2 )( 4 + g 2 u 2 ) 


(D.3) 

(D.4) 


Following the lines of [38] the asymptotic horizontal SU( 2) NLIE variables can 
be determined from the asymptotic Q-functions (1C.ID . Here we just list the final 
formulas: 


b°(u)=b 0 (u-i 7 ), b°(u) = b 0 (u + i 7 ), (D.5) 


where 


2 (g 2 u 2 — 3 i) (1 + 2 i g u + g 2 u 2 ) 
(g 2 u 2 + i) ( g 2 u 2 — 2 i ) ( g 2 u 2 + 3«) ’ 

2 (^ 7 2 u 2 + 3 f) (1 — 2 i g u + g 2 u 2 ) 
(g 2 u 2 — i ) ( g 2 u 2 + 2 i) ( g 2 u 2 — 3 i)' 


(D. 6 ) 

(D.7) 


and 0 < 7 <l/ 2 is the arbitrary contour shift parameter of the horizontal £77(2) 
NLIE. 
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